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Semiconductor heterostructures with prescribed energy de- 
pendence of the transmittance can be designed by combining: 
a) Pade approximant reconstruction of the S-matrix; b) in- 
verse scattering theory for Schrodinger's equation; c) a uni- 
tary transformation which takes into account the variable 
mass effects. The resultant continuous concentration pro- 
file can be digitized into an easily realizable rectangular-wells 
structure. For illustration, we give the specifications of a 2 
narrow band-pass 12 layer Al c Ga\- c As filter with the high 
energy peak more than twice narrower than the other. 
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The inverse scattering method [[1]-|3| allows the recon- 
struction of the potential, V(x), in the one dimensional 
Schrodinger equation, 



2mo 



ip" + V(x)i> = Etj), 



(1) 



from the energy dependence of the S-matrix (complex 
reflection and transmission coefficients). Let us briefly 
review the necessary results, referring to || for details. 
Setting E = h 2 k 2 /2mo > 0, the transmission, T(fc), and 
the left/right reflection coefficients, R±(k) satisfy 

R + (k)T(—k) + R-{—k)T(k) = 0, (2) 
T(k)T(-k) + R±(k)R±(-k) = 1, (3) 
T(k) = [T(-k)]*; R ± (k) = [R ± (-k)}*, (4) 

where * denotes complex conjugation. In addition, 
l-ft±(^)l < 1) except at k = 0, where generically one 
has i?_(0) = —1. Moreover, liim^oo R±(k) = 0. 

The scattering coefficients T(k) and R±(k) can be con- 
tinued analytically to the upper half-plane of the com- 
plex variable k. We will consider only the case when 
V(x) — for x < and there are no bound states. Then, 
for Im(fc) > 0, T(k) and the left-reflection coefficient, 
R—(k), are analytic functions having no zeros. The right- 
reflection coefficient, i? + (fc), may have additional poles. 
The relations (§,[|) remain valid in the whole upper com- 
plex half-plane, while eqs. (Q) hold with — k replaced by 
— k* in the right hand sides. 

The potential V(x) is reconstructed as 



h 2 d 



V(x) = — —K-(x,x), 
mo ax 



(5) 



where the function K_(x,y) satisfies Marchenko's inte- 
gral equation 



K-(x,y) + M-(x + y) + / Af_(y + s)K-(x, s)ds = 0. 

J-v 

(6) 

Here, M-(u) is the Fourier transform of the (complex) 
left-reflection coefficient, R_(k), which is therefore the 
input information required for the inverse technique. 

Although inverse scattering methods have been exten- 
sively used in various branches of physics and engineer- 
ing , they have not yet been used in heterostructure 
design. In this Letter we give a brief outline of the appli- 
cation of these techniques to the design of electron-wave 
filters — semiconductor heterostructures having a pre- 
scribed transmittance, |T(fc)| 2 . So far, such filters had 
been designed using methods adapted from optical filter 
synthesis [pi. These structures may be useful for far- 
infrared lasers and detectors 

In the effective mass - envelope function approxima- 
tion, the electron states in the conduction band of a het- 
erostructure having a position dependent composition (x 
is the deposition direction and c(x) the alloy concentra- 
tion) are solutions of the BenDaniel and Duke equation 



Tr 



2m(x) 



VV> + £ c (x)i> = Etp, 



(7) 



where m(x) = m e //[c(x)] and £ c (x) = S con d[c(x)] are 
the local (isotropic) effective mass and conduction band 
offset of conduction band electrons. We consider here 
only the case when the device is embedded (at least in an 
asymptotic sense) in a crystal with uniform composition, 
i.e. we assume that c = for x < and in the limit 
for x — ► +oo; let be the conduction-band effective 
mass in the embedding material. In many systems of 
interest, such as Al c Ga\- c As, the dependence of m e // 
and £ C ond on c is approximately linear [fill : m e / / (c) = 
a + be, £ CO nd{c) = a + (3c. The mass - concentration 
dependence can then be inverted, to get: 



E c (m) = A + Bin. 



(8) 



For Al c Ga\^ c As we use the numerical values A = 
— 624.1meV and B = 9315meV/m e , where m e is the free 
electron mass. 

After separating away the variables associated with 
motions orthogonal to the x direction in eq. (|7]) and 
setting the transverse quasimomentum qj^ = 0, the en- 
velope function for the x direction, /, satisfies: 
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fi 2 k 2 

+ B[m(x) - moo] / = -^—f. 
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In the above, we have used (0) and set 



fik = T l /2m co E, 



(9) 



(10) 



the energy E being measured from the conduction band 
edge of the embedding crystal. An additional term 
Et((l±) ([rrioo / m(x) — 1] f(x), with E in ([l(]) replaced by 
E - £ t (q_L), where E t (a ± ) = (l/2)h 2 q 2 L /2m 00 , appears 
in the left hand side of (Q) for ^= 0. In Al c Ga\_ c As the 
effect of this corrections is quite small for .E t (qj_) ~ ksT, 
up to room temperature. Eq. (^|) is similar to the 
one-dimensional Schrodinger equation (Q), but with a 
position-dependent mass. This variability has nontrivial 
dynamical effects jl2|]. Noticing that only m(x) enters 
@, our task is now to determine mass (and hence con- 
centration c(x)) profiles giving rise to the pre-assigned 
transmissivity |T(fc)| 2 . Let us highlight the main prob- 
lems and our approach to solving them. 

1. Eq. (^) is a low energy approximation, valid only 
for the sufficiently small energies at which the conduc- 
tion band energy is parabolic. Since we assume to be 
interested only in this range of energies, the transmit- 
tance behavior can be arbitrarily continued to higher 
energies, respecting the analyticity properties of the S- 
matrix. The design problem now has an infinity of solu- 
tions. We will not explore all of them, but we will use 
Pade approximants [l3|| to fit the transmittance speci- 
fications within given tolerances in the required energy 
range. The non-unicity of the approximants will give free 
parameters which can be used in successive developments 
to optimize properties of the device other than the energy 
dependence of |T| 2 . 

2. The transmittance input data we use give only 
|T(fc)| 2 on the real axis, while the solution of the inverse 
problem for (||) requires also the phases of R± . Generally 
speaking, the necessary phases can be reconstructed by 
using dispersion relations as in particle theory (m). Yet, 
since we use Pade (rational) approximants for fitting the 
energy dependence of ITI 2 , we can perform the analytic 
continuation in a simpler way. In fact, in the case when 
no bound states are present, we let 



n(E) = 



(ii) 



be the Pade approximant fitted to the reflectance. V P (E) 
and Q q {E) are polynomials of degree p and q > p+2 with 
real coefficients. The condition q > p + 2 guarantees that 
71(E) = 0(E~ 2 ) for large E and, therefore, the poten- 
tial is non-singular at x = 0. In accordance with the 
analytic properties of the reflectance on the real axis, (Q- 
@), < 71(E) < 1 for < E < oo and 11(0) = 1. The 
latter implies that V p (0) = Q q (0) = 1. From the former 



follows that if V P (E) has real positive zeros (which corre- 
spond to resonances, where the transmittance attains its 
maximal allowed value, unity) their multiplicities have 
to be even. From the structure of eq. @ , we see that it 
is convenient to work in k space, and to compute (jp, q) 
Pade approximants to the left-reflection coefficient 



R-(k) 



P P (k) 
Q q (kY 



(12) 



where the k — values of the polynomials are taken as 
P p (0) = -Q q (0) = -1. This guarantees that i?_(0) = 
— 1. Let k a (respectively, kj) be the zeros of P (Q). To 
satisfy the analytic properties of R-(k) and the condi- 
tion (|J), the zeros of P and Q must have non-positive 
imaginary parts (negative for P). If k a (kj) is a root 
of P(k) = (Q(k) = 0), then -k* (-k*) must also 
be a zero of the same polynomial: zeros of P and Q 
which are not purely imaginary come in pairs. Using 
the root product expansion of the polynomials we may 

write R_(k) = -nLi(l-£)/n^i(l-]|)- For 
real positive energies the reflectance is given by 



K(E) = R-(k)R-(-k), 



(13) 



where we used (^). Thus, taking into account (|l0|), we 
see that the zeros and poles of the analytic continuation 
to the complex E plane of the Pade approximant ([ll]) for 
the real axis reflectance are connected to the zeros and 
poles of the Pade approxima nt (|l2"|) f or the (complex) re- 
flection coefficient by k = ±y / 2m 00 E/Ti. The sign assign- 
ments are unambiguous in the absence of bound states, 
when R- (k) is analytic and has no zeros in the upper half 
plane. The phase reconstruction for the analytic continu- 
ation of the reflection coefficient R-(k) from the absolute 
values on the real axis is therefore achieved. 

Similarly, we seek the transmission coefficient, as a 
(q, q) Pade approximant having the same poles as i?_(fc): 



T(k) 



fc_L <j _ 1 (fc) 
Q q (k) 



(14) 



Here we used T(0) = 0. The polynomial L q ^i(k) is deter- 
mined by its roots and the condition limfe_ >0O T(fc) = 1. 
From (|l^) and (||4|), the transmittance on the real axis 
is given by T(E) = 1 - R_(k)R_(-k) and Q leads to 
the following equation for the zeros of L q _i(k) 

k 2 Lq^(k)Lq^(-k) = Qq(E) - Vp(E) = 0. (15) 

The zeros of L ? _x(fc) are then given via ( JTo| ) (taken with 
non positive imaginary part) by the zeros of the analytic 
continuation in energy of the absolute value of T(E). 
Finally, R+ (k) can be obtained from (^|) and we have fully 
reconstructed the S-matrix from the Pade approximation 
for the energy dependence of the transmittance, in the 
absence of bound states. 
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3. Although similar, (|jj) differs from (|]) by having a 
variable unknown effective mass. 

In Al c Gai- c As, m e ff increases by about 50% from 
the GaAs value in the direct gap concentration range, 
< c < 0.45. There are two alternative ways to tackle 
this problem. One may reformulate the inverse scattering 
method to deal with equation (||). We prefer to follow 
an alternative, two-step approach. First, we solve an 
auxiliary problem: we find a potential for the constant 
mass Schrodinger equation ([!]) such that the correspond- 
ing iS-matrix coincides with the one reconstructed from 
the transmittance, as in point 2 above. In particular, if 
the scattering data are expressed by Pade approximants, 
then the solution of the inverse problem is much sim- 
plified, and it can be accomplished by solving a linear 
system of equations jl|-[|. Second, we will use a uni- 
tary transformation linking constant and variable mass 
Hamiltonians []l5| and derive a differential equation for 
the mass (concentration) profile in (||). 

In fact, let X(£) be a twice differentiable function, sat- 
isfying X(±oo) = ±oo and having a positive derivative 
= A(£) > 0. Here and in the following we use the 
dot to denote derivation with respect to £. Then, let us 
consider the coordinate transformation 



*(0> 



(16) 



which maps the interval (— oo, +oo) onto itself. As shown 
in pjj ], each such change of variables corresponds to a 
unitary operator, U x , in the space of square integrable 
functions, such that in the new coordinate £ the equation 
(III) transforms to 



d fi 



d£ 2m(£) d£ 



f(0 + Vx(0f(0=Ef(0, (17) 



where 



m(£) = moX 2 



v x (0 = v[x(0] 



8m X 4 



(18) 



2XX-5X 2 ), (19) 



and where we omitted for simplicity the £ dependence of 
the derivatives of X. Eqs. (|l^ - ^9|) and (g) coincide if 



V X (S) = B [m(0 - moo] 



(20) 



where m(£) and Vx*(£) are given by ( |18| , |l9j ). This is a 
nonlinear differential equation of third order for X. Not- 
ing that it does not depend explicitly on £, and introduc- 
ing a new dependent variable S(X) — X yields 



2S"S- 3S" 2 = 



8m S 2 



[Bm (S 2 - Si) - V{X)} , (21) 



To be physically feasible, the solution of ( |2l| ) must be 
regular and bounded by 



m < moS' 2 (£) < m, 



(22) 



where m and m are the minimal and maximal achiev- 
able values of the effective mass, and tend to a/ moo/mo 
in the limits X — > ±oo. A necessary condition for the 
latter condition to hold is B > 0, which is validated for 
GaAlAs. The asymptotic value of the effective mass for 
the embedding crystal is a free parameter. Its value must 
be chosen in a range leading to a feasible concentration 
range, satisfying (|22j). 

In practice, we truncate the auxiliary (constant mass) 
potential V(X) to a finite interval (0, L). For zero po- 
tential Eq. ( pi| ) can be explicitly solved: S(X) = 
[a + bc KX + ce~ KX ] ~~ , where fi 2 n 2 = SBm^ and the 
constants a, 6, c satisfy a 2 ~ Abe = mo/moo. Now, let 
S(X) be a regular solution on one of the intervals (-co, 0) 
or (L,+oo). It satisfies 



S'(X) = ±kS(X) 1 - y/mo/m^SiX) 



(23) 



The prime denotes derivation with respect to X. 



at all points of the respective interval, with the minus 
sign on the former and the plus on the latter. Thus, 
in order to be continuous and achievable, a solution of 
( pl| ) on (0, L) must satisfy (|23| ) at the ends of the in- 
terval. This is a nonlinear Sturm-Liouville problem that 
can be solved by the shooting method. The new coordi- 
nate can be obtained then by integrating A(£) = S(X). 
Reasonable precision may also be obtaineded by using a 
quasi-classical approximate solution to j2l|). 

Thus, we obtain a continuous distribution of mass 
(concentration) yielding the desired transmittance. Tech- 
nological constraints make desirable concentration distri- 
butions in the form of a sequence of steps of digitized 
length (in lattice constants) and height (linear combina- 
tions of several concentrations). This can be achieved 
by approximating the continuous concentration profile 
and by using standard optimization techniques. The 
brute force approach to the inverse problem is to find 
the transfer matrix of (^|) for concentration profiles made 
of discrete steps (parameterized by concentrations and 
lengths) and constrain the resultant transmittance to the 
predetermined profile. This suffers from the fact that the 
corresponding optimization problem is non-convex and, 
even using the largest computers, satisfactory results are 
not guaranteed. To the contrary, our approach offers 
a good starting point for the optimization process. A 
modest exploration of its neighborhood usually provides 
superb results. 

As an illustration, we give in Table | the data for a 
12 layer Al c Ga\- c As filter with digitized steps in width 
and concentration. Quite remarkably, this structure ac- 
complishes the goal of having a lower energy resonance 
more than twice wider than the high-energy one. In Fig. 
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|l| the transmittance computed from (|^) is compared to 
a constant effective mass approximation using (Q) with 
too = Woo. The latter's errors become important at the 
higher resonance, whose width is overestimated by a fac- 
tor of two, thereby losing the main feature of this struc- 
ture! The computation of the continuous profile takes 
several seconds on an IBM RISC 6000 workstation. The 
digitizatization and final optimization were done manu- 
ally. 

Spatial variations of the band gap have been crucial 
for the success of man-made heterostructures. We have 
outlined a technique which takes into account also the 
spatial variations of effective mass, and can be employed 
in the design of various electronic devices. 
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TABLE I. 12 layer 


digitized Al c Ga\- 


c As filter. 


Layer 


Width 


Width 


Al 


# 


(unit a ) 


H 


concentration b 


1 


6 


1.696 


C1+C2 + C3 


2 


9 


2.543 





3 


18 


5.088 


C3 


4 


5 


1.413 


C2 


5 


10 


2.827 


Cl + c 2 


6 


10 


2.827 


Cl 


7 


14 


3.957 


Cl + c 2 


8 


14 


3.957 


Cl 


9 


13 


3.675 


Cl + c 2 


10 


13 


3.675 


Cl 


11 


14 


3.957 


Cl + C 2 


12 


11 


3.109 


cl 


a monoatomic layer width 


= 0.2827 nm 




b ci = 


0.05714, c 2 = 2ci, 
1 , — , * , 


cs = 4ci, bulk concentration = c 2 
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FIG. 1. Energy dependence of the transmittance of 12 layer 
digitized Al c Gai~ c As filter: continuous line — Eq.([|); dotted 
line — constant mass approximation Eq.(|l]) with mo = m^. 
Insert: potential energy profile. 
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